Determination of the tunnel rates through a few-electron quantum dot 
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We demonstrate how rate equations can be employed to find analytical expressions for the 
sequential tunneling current through a quantum dot as a function of the tunnel rates, for 
an arbitrary number of states involved. We apply this method at the one-to-two electron 
transition where the electron states are known exactly. By comparing the obtained expressions 
to experimental data, the tunnel rates for six transitions are extracted. We find that these 
rates depend strongly on the spin and orbital states involved in the tunnel process. 



1 Introduction 

Quantum dots in semiconductor heterostructures have proven to be excellent systems to study 
the charge and spin degree of freedom of single electrons as well as coherent many-body 
processes like the Kondo effect*^. Quantum dots that are defined electrostatically in a two- 
dimensional electron gas are especially interesting, since here both the electron number and the 
tunnel coupling to the leads are fully tunable. The experiments in these systems are currently 
at the level of studying and manipulating the coherence properties of single electrons. Besides 
being interesting from a fundamental point of view, coherent control is an essential element of 

Q 

proposals to use the spin of a single electron as a building block for a quantum computer . 

In order to benefit fully from the tunability of the tunnel barriers, one needs to know the 
tunnel rate for all possible charge transitions. It was recently demonstrated that the tunnel 
rate through a single barrier can be determined by applying voltage pulses to one of the surface 
gates and monitoring the response of the quantum dot with a charge detector!^. This method 
can also give the rates of transitions involving excited states, but, due to the bandwidth of 
most measurement setups, its applicability lies mainly in the regime of very weak coupling to 
the leads (tunnel rate < 100 kHz). For stronger coupling, the tunnel rates can in principle be 
determined from the measured current as a function of the source-drain bias voltage, by solving 
rate equations. This has been demonstrated for a system with tw o allo wed transitions, both in 
the limit of very asymmetric and of completely symmetric barriers^^^^. 

In this work, we demonstrate that an expression for the current through the dot as a function 
of the tunnel rates can be found for an arbitrary number of possible transitions using rate equa- 
tions. Then, we apply this method at the l<->2 electron transition and give exact expressions 
for the case of symmetric barriers. By comparing these expressions to measurements, we deter- 
mine the tunnel rates of the six lowest-energy transitions. Finally, we comment on the observed 
dependence of the tunnel rate on the orbital and spin state involved in the tunnel process. 



2 Calculation of the current using rate equations 

2.1 Model and assumptions 

The system we study is a quantum dot weakly coupled to a source and a drain reservoir. We 
calculate the sequential tunneling current flowing through the dot, thus neglecting higher-order 
tunneling events. Also, we assume that relaxation processes within the quantum dot are much 
slower than the tunnel rates. This is generally true for relaxation between states with different 



spin (which w e wi h consider in this work), but not necessarily for relaxation between states with 
the same spin^^El, 

Since in our experiments both the energy spacing of the quantum dot states and the source- 
drain bias are much larger than the electron temperature, we take T = 0. Thus, all states below 
the electrochemical potentials of the reservoirs are filled and all states above are empty. This 
implies that when a bias voltage is applied between the source and the drain, the electrons can 
only tunnel onto the dot from the source contact and can only tunnel off the dot to the drain, 
i.e. they tunnel in the forward direction only (exceptions will be explicitly stated). 

We focus on the case where the tunnel rates through the source and the drain barrier are 
equal, since this simplifies the resulting equations considerably. In the experiments, the barriers 
can easily be tuned to make this assumption valid. 



2.2 Calculation of steady-state occupation probabilities 



To find the current flowing through the dot, we first calculate the occupation probabilities for 
all relevant states by solving the rate equations for the available transitions. Consider a dot 
where energy states contribute to charge transport. The time evolution of the occupation 
probability po of an energy state |0) depends on the occupation probabilities of the other energy 
states and the tunnel rates between the states: 
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where Tj^i is defined as the tunnel rate between the initial state \i) and final state |/). We can 
write down the rate equations for all N+1 levels conveniently within a matrix representation: 
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We call r the transition matrix. The first column of this matrix (except for the diagonal element) 
consists of the probabilities per unit time, T^fi, that the dot makes a transition to state \k) if 
it was initially in |0). These probabilities follow Poissonian statistics. The next N columns 
correspond to the dot initially being in . . ., |A^), respectively. The diagonal elements of the 
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matrix, Ta depend on the non-diagonal elements via the relation, Fj 

We are interested in the steady-state solution of Eq. [2 Since in this case the occupation 
probability distribution is constant in time, we find the steady-state occupation probability 
distribution by solving 
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i.e. determine the eigenvector of T with eigenvalue 0. 



2.3 Calculation of the current from the occupation probabilities 

The current can be calculated from the steady-state occupation probabilities by multiplying 
these probabilities with the tunnel rates through one of the barriers. The currents through 
the two barriers are equal in the steady state. Consider transitions between states with 
M electrons on the dot, having occupation probabilities Pi(M)^ states Ijm+i) with M+1 
electrons on the dot, having occupation probabilities Pj(M+i)- If only forward tunneling is 



possible, the current through one barrier wih be due to electrons added to the dot when it 
holds M electrons, whereas the current through the other barrier is carried by electrons that are 
extracted from the dot when it is in one of the (M+l)-electron states. The current is therefore 



^ = PiW rj{M+l),i(M) = e XI Pj{M+l)Y.'^iW,jiM+l) (4) 
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If electrons can also tunnel in the opposite direction (reverse tunneling), the current can still be 
calculated in a similar way by making suitable substitutions, as shown below. 



3 Expressions for the current at the 1^2 electron transition 

3.1 Energies and electrochemical potentials of the 1<^2 electron transition 

We apply the method explained in Section 2 to the transition between 1 and 2 electrons on 
the dot and find expressions for the current for different source-drain bias conditions. First we 
write down the energies of the relevant states and from these we deduce the electrochemical 
potentials of all the allowed transitions. To observe spin effects, we assume that a large in-plane 
magnetic field B// is applied to the dot, which has a negligible effect on the orbitals but causes a 

large Zeeman energy splitting AEz= gfisB/^ between the spin states^. The two lowest-energy 
one-electron states are the two spin states of the lowest orbital, ||) and ||), with energies 
and Ei = E^+AEz respectively. 

The ground state for two electrons on the dot is a spin singlet (total spin quantum number 
5 = 0)1^, formed by the two electrons occupying the lowest orbital with their spins anti-parallel: 
I -S*) = (I Ti ) — I iT))/'v/2- The first excited states are the spin triplets (5=1), where the antisym- 
metry of the two-electron wave function requires one electron to occupy a higher orbital. The 
three triplet states are degenerate at zero magnetic field, but acquire different Zeeman energy 
shifts in finite magnetic fields because their spin z-components (quantum number ms) differ: 
ms = +l for |r+) = |TT), ms = for jTo) = (| n) + UT))/\/2 and ms = -l for |r_) = |ii). We 
can write the energies of the two-electron states in terms of the single-particle energies of the 
two electrons plus a charging energy Ec which accounts for the Coulomb interactions: 

Es = E-^ + E^ + Ec = 2E^ + /S.Ez + Ec 
Etj^ = 2E^ + Est +-E'C 

Eto =Ei^ + Ei+EsT+Ec = 2E^+EsT + AEz+Ec 
Ej'_ =2Ej^-\- Egx~^ Ec = ^E^ -\- Egx^'^^Ez ~\~ Ec , 

with Est denoting the singlet-triplet energy difference in the absence of Zeeman splitting. 

The electrochemical potentials fi for the different transitions are now easily calculated using 
f^a^b = El, — Ea This yields for the l<->2 electron transitions: 

E^ + Ec + AEz (5) 

E^ + Ec (6) 

E-^ + Ec + Est (7) 

Ei^ + Ec + AEz + Est- (8) 

We have omitted the transitions |^T_ and J,<->-T+, since these transitions require a change in 
the spin z-component of more than | and thus a spin-fiip is needed in the tunneling process. 
Since the single-spin Zeeman relaxation'^ as well as the triplet-to-singlet relaxation!^ is very 
slow, we can neglect these spin-blocked transitions. 
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3.2 Solving the rate equations of the 1^2 electron transition 



By adjusting the electrochemical potentials of the source, ixs, and the drain, ^u/j, relative to the 
ladder of electrochemical potentials in the dot, six different configurations can be obtained where 
current flows. These configurations are depicted in Figs, ^-f, and we label them A through 
F, respectively. (Note that current can only flow if the transition between the one-electron and 
two-electron ground states, t'^'S', is in the bias window, i.e. > fJ'i^s > A'D-) 
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Figure 1; (a)-(f) Electrochemical potential diagrams at the 1^2 electron transition, corresponding to the con- 
figurations A-F respectively. Some transitions (indicated by arrows) are possible even when the corresponding 

electrochemical potential is not in the bias window. 



To obtain expressions for the current, we first identify which transitions are possible for a 
specific configuration and set the tunnel rate for all other transitions to zero. Also, we take into 
account that for some configurations electrons can tunnel in both forward and reverse direction, 
which yields a doubling of the transition probability. Then, we find the steady-state occupation 
probabilities of all the states in terms of the tunnel rates by calculating the eigenvector of the 
resulting transition matrix with eigenvalue zero. This is done using the program Mathematica. 
From the occupation probabilities we obtain an analytical expression for the current as a function 
of the tunnel rates, as explained before. We now demonstrate the construction of the transition 
matrix for configurations C and F. 

In the basis p = {pi,Pi,ps,PT+,PTo,PT-), the full transition matrix for the 1^2 electron 
transition is 
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In configuration F, all transitions are possible and electrons can only flow in one direction (see 
Fig. W^), and therefore Eq. IHl directly gives us the relevant matrix. 

In configuration C, the transition from \S) to ||) is blocked, and the transition from ||) 
to IS") is possible via both barriers. Thus, we get the transition matrix for C by making the 
substitutions F|5 and Tsi 2F51 in Eq. |^ 

The transition matrices for the other configurations can be derived in a similar way. The 
calculated expressions for the current in configurations A through F are given in the Appendix. 



4 Extraction of the tunnel rates from experimental data 

In this section, we use the expressions for the current obtained in the previous section (and given 
in the Appendix) to extract the tunnel rates for all six transitions from experimental data at B n 
= 12 T. This data was previouslv analyzed in the context of spin filtering. Experimental details 
can be found in Ref.'^. In Figs. |^ andUJa we present traces of the current through the dot as 
a function of gate voltage Vq- Changing Yq shifts the whole ladder of electrochemical potentials 
and thus allows us to scan through different configurations as indicated by letters A through F . 
The two traces in Figs. 0^ and|l]3 are taken at different source-drain bias voltage Vsd- Due to 
capacitive coupling of the reservoirs to the dot, changing Vsd also slightly changes the barriers 
and thereby the absolute values of the tunnel rates. However, we assume the relative values of 
the tunnel rates to be independent of Vsd and thus the ratio between the rates to be constant. 




Figure 2: (a)-(b) Current / as a function of Vb at B// = 12 T for (a) Vsd = 310 /iV and (b) Vsd = 800 /xV. 
In these graphs, A-F indicate the configuration of the levels, (c) The current / via the singlet state (S) and via 
both the singlet and all triplet states (S+T) as a function of Bx_ (Vso=l meV). The currents are normalized to 
the values at 5^=0.4 T. The Zeeman energy is not resolved at these low magnetic fields. 

From the trace shown in Fig. we find that the current in configuration A is 0.41 pA, 
yielding the tunnel rate T s\ = 5.1 MHz. The measured current in configuration D is 0.77 pA, 
which leads to Tsi = 9.3 MHz, using F^^ = 5.1 MHz. We note that there is a spin dependence 
in the tunnel rate: « 2F5'^. This spin dependence has been observed before at the 0<->l 
electron transition in the same device^l, and might be attributed to exchange interactions in 
the leads close to the dot. The ratio for spin- up and spin-down electrons tunneling to the same 
orbital, a, is assumed to be constant around the 1 <-> 2 electron transition and independent of 
Vsd and Vg- In the rest of this section we use a = 2. 

The current in region A is 0.44 pA close to region B, and the tunnel rate Ts], is calculated 
to be 5.5 MHz. It differs slightly from the value close to region D obtained above because the 
barriers are not completely independent of Vq- In region B, the measured current is 14 pA, 
yielding F^+t = 0.27 GHz. 

We now turn our attention to Fig. EJj. Here, the current in configuration D is Id = 0.60 pA, 
which gives F^^ = 3.8 MHz, Tsi = 7.6 MHz and F^^i = 0.18 GHz, assuming the ratio of the 
rates is fixed. For the current in configuration E we find Ie = 13 pA. From this number we 
deduce T^q^ = 0.09 GHz and Tx^i = aTT^^ = 0.18 GHz. From the magnitude of the current in 
F, If = 14 pA, we find Ft.; = 0.18 GHz. 

Since all the rates have been determined, we can predict the current in configuration C. The 
rates determined above predict a current of 9.2 pA, in good agreement with the measured value 
of 10 pA. 

We finally comment on the strong dependence of the tunnel rates on the orbital state involved 
in the tunnel process: tunneling to a triplet state is more than an order of magnitude faster than 
tunneling to the singlet state. This effect can be explained by the different spatial distribution of 
the singlet state and the triplet states. In a singlet state, both electrons are in the lowest orbital 
state, and are therefore mainly located near the center of the doiJ^. In a triplet state, one electron 
occupies an orbital excited state which has more weight near the edge of the dot, and therefore 



is expected to have a larger overlap with the reservoirs. To illustrate the significance of the 
spatial weight, results are presented from a measurement on a similar device in a perpendicular 
magnetic field B±. The current via the singlet state and the current via both the singlet and 
the triplets are plotted as a function of B± in Fig. |^. The perpendicular field causes the wave 
functions to shrink in the plane of the two-dimensional electron gas, reducing the overlap with 
the leads and thereby reducing the corresponding tunnel rates. The decrease in current for 
increasing B± is clearly observed. 
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Appendix 



The analytical expressions for the current as a function of the tunnel rates in configurations A 
through F, I a through Ip, are: 

Ia = ^eFsT (10) 
Ib = ^e(rsT + rT+T) (11) 
= / 8 " 2 (12) 



Id = ^e{Ts^+Tsi) (13) 
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